Effective Action Studies of Quantum Hall Spin Textures 
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We report on analytic and numerical studies of spin tex- 
tures in quantum Hall systems using a long-wavelength effec- 
tive action for the magnetic degrees of freedom derived previ- 
ously. The majority of our results concern skyrmions or soli- 
tons of this action. We have constructed approximate analytic 
solutions for skyrmions of arbitrary topological and electric 
charge and derived expressions for their energies and charge 
and spin radii. We describe a combined shooting/relaxational 
technique for numerical determination of the skyrmion pro- 
files and present results that compare favorably with the an- 
alytic treatment as well as with Hartree-Fock studies of these 
objects. In addition, we describe a treatment of textures at 
the edges of quantum Hall systems within this approach and 
provide details not reported previously. 
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I. INTRODUCTION 

Following the demonstration in that quantum Hall 
ferromagnets contain skyrmions in their quasiparticle 
spectra, there has been a great deal of work on spin tex- 
tures in quantum Hall systems. This work has received 
considerable impetus from a set of experiments j|-[D that 
have adduced evidence for skyrmions being the lowest en- 
ergy quasiparticles at v = 1 as well as for their not being 
so at v = 3 and higher odd integer fillings ||[7|- The 
skyrmions themselves have been investigated by several 
authors ||^| and a skyrme crystal phase has been stud- 
ied pc| ] with a recent experiment appearing to find a 
transition into it at a finite temperature pd[ |. In addi- 
tion, textured quasiparticles have been discussed in the 
context of double layer systems and invoked in ex- 
plaining a novel phase transition in the latter [|l3| . More 
recently, it has been shown that ferromagnetic quantum 
Hall states can exhibit textured edges @,^5| ■ 

The analysis in [jjj was based on a long wavelength 
effective action for the spin dynamics, supplemented by 
a charge density-topological density constraint, that al- 
lowed an analytic computation of the properties of the 
skyrmions in the limit of small Zeeman energies. No de- 
tails of this computation were given there and one of the 
purposes of this paper is to provide them. A second is 
to report numerical studies of the effective action, that 
have a region of validity beyond the analytic approxi- 
mations, and compare extremely well with those using 
the Hartree-Fock technique of Fertig et. al. M . We also 



report here on the details of effective action studies of 
the textured edges mentioned in (m). We note that the 
numerical studies have been described previously in a 
thesis by one of us |l6| ] and are being reported here for 
ease of access. In the interim, effective action results on 
skyrmion properties have been reported by Abolfath et. 
al. using methods similar to ours |17j] , and by Rao et. al. 
using a variational approach |Q , while Moon and Mullen 
have presented an improved action that is accurate even 
for small skyrmion sizes Jf^] . 

The paper is organized as follows. After introducing 
the effective action in Sec. ||, we outline the analytic 
derivation of skyrmion properties at small Zeeman ener- 
gies m Sec. |H In Sec. [Fv|, we discuss the technique 
utilized in their numerical computation. In Sec. [v|, we 
present the numerical results for charge I and charge 2 
skyrmions at v = 1, and compare to Hartree-Fock cal- 
culations and analytic expansions. We then discuss the 
asymptotic shape of the skyrmions (Sec. VI). This is 



followed by results for textured edges (Sec. VII) and a 
brief summary (Sec. VIII). 



II. EFFECTIVE ACTION AND SKYRMIONS 

In the long wavelength limit, the spin degree of free- 
dom of a quantum Hall (QH) ferromagnet is described 
by a unit vector n(r). A distinguishing feature of QH 
ferromagnets is that the topological density of the spin 
field, q(r) = n ■ (d x ii x d y n)/4:Tr, is proportional to the 
deviation of the charge density, p(r), from its background 
value, p — p = v FM q where p — v FM l2-Ki 2 . {v fm is the 
filling factor of the ferromagnetic component of the QH 
liquid {e.g. v FM = 1 at v = 3), and I — yJhc/eB is the 
magnetic length.) This enables one to formally integrate 
out the dynamics of the charge and obtain an effective 
Lagrangian for the system in terms of the spin field alone 
whose coefficients can be fixed using known long wave- 
length quantities 0: 



L eff = -pA{n)d t n 



:Ps(Vn) 2 



4* % I drr 



,9( r M r ') 



(1) 



Here, A is the vector potential of a unit magnetic 
monopole, p s is the spin stiffness |2(| and e the dielectric 
constant of the background semiconductor. Technically, 
Ceff describes an 0(3) a— model. 
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As we will be interested in this paper in static config- 
urations of the field, we will seek to minimize the energy 
functional 



E = - d 



i c eff ~ -pA{n)d t n] = E G 



E 2 



Er 



(2) 



where the gradient {Eg), the Zeeman (Ez) and the 
Coulomb (Ec) energies can be read off from Eq. ([!]). 
In the following we will measure all energies in units of 
e 2 /e£ and lengths in units of I. 

For finite energy configurations, the finiteness of the 
gradient term requires that the field approach a common 
value at infinity regardless of direction. In such cases, the 
plane can be compactified to a sphere and n(r) gives a 
map from Si (the compactified plane) to S2 (spin space 
or the target space of the o— model). Such maps are 
characterized by an integer topological charge Z, which 
can be expressed as Z = J d 2 rq(r), where q is the topo- 
logical (Pontryagin) density introduced earlier. As noted 
there, the topological density plays a crucial role in the 
physics of QH ferromagnets as it is proportional to the 
charge density of the underlying itinerant system. Conse- 
quently, localized configurations of topological charge Z, 
which we shall term charge Z skyrmions, carry electric 
charge v FM Z\e\. 

In the absence of the Zeeman and Coulomb terms, the 
energy functional (|J) is scale invariant and rotationally 
invariant but on account of the non-linearity implicit in 
the definition of n, finding its skyrmion solutions is a 
non-trivial problem, solved previously by Belavin and 
Polyakov (BP) The additional terms break both 

symmetries although they preserve rotations about the 
field axis which we shall take to be the z-axis. For the 
full action, the BP technique (based on achieving a Bogo- 
mol'nyi bound |^3| for the action) breaks down with the 
non-locality of the Coulomb term making matters even 
worse. Consequently, it is necessary to attack the Euler- 
Lagrange equation for the energy functional directly. 

We now derive this equation for the skyrmion con- 
figurations by making an ansatz that has topological 
charge Z and depends on one unknown function. Cal- 
culating the energy of this configuration and minimizing 
with respect to the unknown function gives an integro- 
differential equation. 

In the ground state the spin of the ferromagnetic com- 
ponent is polarized: n = z. For a skyrmion with topolog- 
ical charge Z ', we make the ansatz (in polar coordinates 
(r, 9) with their origin at the center of the skyrmion) : 



Anr dr 



(4) 



For a finite energy configuration, the spin vector must be 
aligned with the external magnetic field at infinity. At 
the center the spin points in the opposite direction. This 
leads to the boundary conditions /(0) = —1, /(oo) = 1, 
and (|J) shows that (||) describes configurations with 
topological charge Z. Note that the topological density 
is independent of 8, i.e. it is spatially rotationally in- 
variant; since a shift 9 —> 8 — o in (j^) is a spin rotation 
it is clear that the topological density is invariant under 
spin rotations about the z-axis as well. Note however, 
that the spin field itself is not separately invariant under 
rotations and spin rotations — it is invariant only un- 
der the combination generated by L z + ZS Z . It is easy 
to convince oneself that respecting both symmetries is 
incompatible with non-trivial skyrmion solutions. Con- 
sequently, our ansatz is the maximally symmetric one 
possible. Examples of spin configurations for Z = 1 and 
Z = 2 are sketched in Fig. |l|. 
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n x = v 7 ! -P(r) cos(Z6) , 
n y = y/l-p(r)sm(ZO) , 
n z = f(r) ■ 

This leads to the topological density 



(3) FIG. 1. Skyrmion spin vector field, n(r), and its projection 

onto the plane of the electron gas for Z = 1 (top figure) and 
Z — 2 (bottom figure). 
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The function f(r) in (||) is determined by minimiz- 
ing the energy (^). We now calculate the energy for a 
skyrmion of the form (|^). The gradient and Zeeman en- 
ergies are: 



2 1-/ 2 , rf 2 



E G = 7rp s / dr[Z 
Jo 

1 f 00 
Ez = --^v PM g J drrf(r) , 



r I-P 



(5) 
(6) 



where g = gpBB/(e 2 /e£) and /' = df /dr. We write the 
Coulomb energy as 



E c = -v 2 FM Z 2 ^ drfVc(rJ') , 



(7) 



where 
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Vc(r, f) = -^J dr'f(r') J dO'[r 2 + r' 2 
-2rr' cos(0' - 9)]- 1/2 



r + r 2 r + r 



(8) 



Here, F(tp,x) is the elliptic integral of the first kind P4| . 

Adding all energy contributions, (||), (0) and (0), we 
obtain the final expression for the total energy: 



E 



dr[irp s (Z 



2 W 2 



l-f 



-) 



\v FM ~grf + \v 2 FM Z 2 fV c (r,f)\ 



(9) 



Varying the energy functional with respect to / we find 
(observing that SEq/SJ gives a contribution —2dVc/dr) 
the non-linear and non-local integro-differential equation: 



2 f" + rf}(l-f 2 )+r 2 ff 2 + lZ 2 f + 



gv F M 2 

r 



v 2 Z 2 dV c 



4:-Kp s dr 



(i-/T = o 



(10) 



Note that Vc is a functional of /. 

As we noted earlier, the gradient energy, Eq, is scale 
invariant and that minimizing it alone leads to the BP 
solutions of arbitrary size A. Rewriting Eq. (||) as 



E G = ^Anp s Z + np s 



i-/ 2 



we find the Bogomol'nyi bound on the energy 
E G [f] >4*P,\Z\ , 



(11) 



(12) 



and solving rf ± Z(l — f 2 ) = we obtain the BP solu- 
tions, 



m = 



(r/X) 



2\Z\ 



(13) 



with energy Eq = 4%p s \Z\. 

The energy E is independent of the sign of Z, hence 
each solution /(r) to ( [To| ) describes both a skyrmion (Z > 
0) and an antiskyrmion (Z < 0). In the following, when 
we talk about skyrmions, and assume Z is positive, the 
results are equally valid for antiskyrmions with negative 
charge provided we replace Z by \Z\. 

By contrast to the BP solutions, the skyrmions in the 
QH fcrromagnet have a definite size set by a competi- 
tion between the Coulomb energy, Eq, and the Zeeman 
energy, Ez ■ The Coulomb energy favors large skyrmions 
since it decreases when charge is spread out. The Zeeman 
energy, on the other hand, increases when more spins are 
flipped, i.e. when the skyrmion becomes larger, and thus 
favors small skyrmions. This leads to a competition that 
determines an optimal size and energy of the skyrmion. 

In the next section we will construct approximate an- 
alytic solutions to Eq. ( |l0| ) and in sections IV and ^ we 
will describe numerical solutions obtained by a relaxation 
procedure. 



III. ANALYTIC TREATMENT 

We begin by recalling the Bogomol'nyi bound (|l^) on 
the gradient energy of a configuration with topological 
charge Z, Ec[f] > Anp s Z. The Coulomb self-interaction 
of a charge distribution is non-negative, Eq > 0, and 
upon dropping a negative extensive constant which im- 
plies that we are measuring excitation energies relative 
to the ferromagnetic ground state configuration, Ez > 
as well. Consequently, E[f] > 4:irp s Z. 

We will be interested in large skyrmions, for which the 
long-wavelength effective action approach will be accu- 
rate. As remarked previously, this will be the case when 
Eq "dominates" Ez, i.e. as g — » 0. In this limit of diver- 
gent skyrmion size, it is clear that Eq <C 1 whence by a 
general balance argument the value of Ez ~ Ec < 1 as 
well. Consequently, as g — > 0, the energy of the skyrmion 
will approach the Bogomol'nyi bound and we expect that 
the solution itself will converge to a BP solution with an 
appropriately chosen A. 

It turns out that this observation is sufficient to deter- 
mine the leading small g characteristics of the skyrmions 
with Z > 2, but needs to be supplemented by global con- 
siderations for Z = 1. These follow from identifying the 
various length scales in the problem by pairwise balanc- 
ing the terms, Eq ~ 1, Eq ~ 1/Aandi?z ~ ,gA 2 , in the 
energy functional. This yields three scales: Ri ~ 1 (Eq 
and E c ), i? 2 - l/g 1/3 (Ez and E c ) and R 3 - l/g 1 ' 2 
(Ez and Eq). Note that two of these scales diverge as 
g — > with R 3 ^> R 2 . Hence a global solution to (12) in 
this limit can be expected to exhibit both of these scales. 
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We have not succeeded in constructing such a global 
solution, which appears to be a difficult task due to the 
non-locality of the equation. Instead we will attempt 
to formally perturb around the BP solutions. To lowest 
order, this is a problem in degenerate perturbation theory 
and requires merely that we minimize Eq + Ez + Ec in 
the subspace of the BP solutions. By construction, Eq 
is a constant in this subspace. The relevant integral for 
the Zceman term can be done analytically, and yields, 



E 7 . = v fm -k 



2 2 / z - 1 csc(n/Z) _ 



Z 



(14) 



For Z — 1, the double integral for the Coulomb term can 
be carried out as well: 



E r = vi 



3tt 2 1 



" " FM 128A 



(Z = l) 



(15) 



However, for Z > 2 we had to take recourse to numerical 
integration with the results, 



E c = v_ 



1.49 

FM ^ 

4.16 



FM 

„ 8.41 



A 



(Z = 2) 
(Z = 3) 
(Z = A) . 



(16) 



Note that Ez (Z = 1) is infinite. Putting this aside 
for the moment, we minimize Ez + Ec with respect to A 
for Z > 2 to find, 

A = Q.780(v FM /g) 1/3 and 

E = 2 [^+ 1 M4m~9) 1/3 } {Z = 2) 
A = l.29(v FM /g) 1/3 and 

^ = 3[^+1.61(^ M 9) 1/3 ] (Z = 3) 
A = \.7h(v FM /~g) 1/3 and 

E = 4[^+1.80(4 M g) 1/3 } (Z = 4) . (17) 



We will see later on that these provide an excellent de- 
scription of the numerically determined solutions at small 
g, as expected. 

The physically interesting case of Z = 1 is more diffi- 
cult. One can see from Eq. ( [l3|) that 1 — f(r) ~ 1/r at 
large r, whence Ez diverges logarithmically with system 
size. Consequently, the BP solutions cannot be used to 
globally approximate the true solutions at any value of 
g ^ 0; the latter must exhibit a faster decay of 1 — f(r) 
at large r. 

The asymptotic behavior of the true solutions can be 
obtained directly from the equation of motion (jl^) . At 
points far from the core, the Zeeman term in the equa- 
tion dominates over the Coulomb term. The latter is 



the potential of a localized charge distribution and will 
decay asymptotically as 1/r. Dropping the latter, set- 
ting / = 1 — ip 2 /2 and linearizing shows that ip satisfies 
Bessel's equation, 



+ r ^> _ ( Z 2 + ^^-U = 

47T/9 S 

with the physical solution (ip — > as r - 
Kz(r/Rz) in terms of the scale Rz - 



oo 



(18) 

ip cx 
introduced 

VFM9 

earlier. As Kz(x) ~ \p^ e ~ x f° r x ^S> 1, the true asymp- 
totic behavior of 1 — / is exponential for all Z, different 
from the corresponding BP solutions. 

For Z > 2, this discrepancy will (presumably) be at- 
tenuated by going to higher orders in the expansion about 
the BP solutions - but it has no bearing on the validity 
of the leading small g behavior derived previously. For 
Z = 1, we will make the assumption that matching an 
inner BP solution to an outer exponentially decaying so- 
lution (the minimal fix for the divergence encountered) 
will yield the leading small g behavior. (Such a solution 
will certainly converge pointwisc to a BP solution, but 
we have not shown that the corrections to it are sub- 
dominant. For Z > 2 one can readily show that this is 
the case by linearizing about the optimized BP solution 
p5[.) More precisely, we assume that, 



f(r) 



(r/A) 2 



(r/X) 2 

-,2 



r < R 3 



1 



2 Kli R 3 } 



r > A 



(19) 



is the desired solution at small g (for Z = 1: a\ = 
16A 2 / . It is easy to check that with this choice of 
a, both pieces of the definition of / agree on the interval 
A <C r <§; i?3 where they overlap. 

Assuming A <^ i?3, we find that, with the above choice 
of/, 



Ez =4v FM g\ 2 ln(^-) 



(20) 



where ln(ai) = | + 7, i.e. ai = 2.9365. Combining this 
with our earlier expression for Ec, we are required to 
minimize 



(21) 



with a = 4v FMl (3 = ^ 1 ^ M and S = 47rp s / 'a\v FM in 
the region A <§; R3, i.e. X^/g <C 1. Differentiating with 
respect to A and setting A 3 — (3/3/2a)A' 3 then leads to 
the transcendental equation, 



1 

A 73 " 



(22) 
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with c = (2a/3/3)(<5/e) 3 / 2 . As c involves p s , which de- 
pends on v, it does not scale trivially with v\ for v — 1 
it equals 0.0178. Eq. (22) has solutions only for g 



below some critical value, g C3 . For v = 1, one finds 
g C3 = 4.3 x 10~ 5 . At this point A = 12.6, hence the 
condition \\fg < 1 is obeyed and we also note that 
c/^= 2.7. 

Equation (E2T) can be solved perturbatively at small g 
by replacing A^ in the argument of the logarithm by the 
full RHS. This yields, to second order in this procedure, 



1 



2a 



g {ln(c/ y^) + hiln(c/ y^) 
0(lnln(c/V5)/ln(c/v^))} 



A 3 3/3 



for the parameter A and thereafter 



(23) 



0((lnln( C /v^)/ln( C /Vg)) 2 )} (24) 



2\n(c/Vg) 

(with A = 3 * (a/3) 1 / 3 (/3/2) 2 / 3 = 0.7838) for the energy 
of the skyrmion. In all of these expressions, the scale for 
g is set by c 2 « 3 * 10~ 4 which is quite small. This comes 
about as the argument of the logarithm is essentially the 
ratio R3/R2 and the asymptotic regime starts to make 
sense only when this ratio is O(l). As this ratio grows 
extremely weakly, as g 1 ^ 6 , one has to go to fairly small 
g before entering asymptopia and as a result the scale of 
the logarithms, which is set by the scale of this crossover, 
is small. We also note, as discussed at more length else- 
where, that the logarithmic enhancement of the energies 
of the Z = 1 skyrmions relative to their Z — 2 cousins 
causes them to bind at extremely small values of g 

Finally, we can use the leading solutions developed 
here to compute various measures of the size of the 
skyrmions. Quantities of particular interest are the spin 
of the skyrmions and their root mean squared spin and 
charge radii defined via: 

1 



Sz = pJ <f r-(l - /(r)) 
R 2 s = {r 2 ) s = p [<frr 2 l(l-f(r))/S z 



R. 



(r 2 



d rr q(r) , 



(25) 



where we have approximated the local density of elec- 
trons, p(r) = p + u FM q(r), by p = v fm /2tt. These 
quantities are not all independent; an integration by 
parts shows that S z = v FM R 2 /2, and, moreover, since 
Ez = gS z , the spin is obtained from ( |l4] ) and (pc|). 
For the spin radius we find, 



R, 



R3 



3 kn 1/2 (A) 

v ai A * 



(Z = l) 



=2 \/^ Ainl/2( 5 ) {z=2) 



IT 



sec 



(V^)) 1/2 A (Z>3) 



(26) 



Note the feature, that de- 



with ln(a 2 ) = 7 - ± - 
pending on the charge of the skyrmion and the highest 
power of radius in the integrand, we get different scales 
entering the leading dependence. 



IV. NUMERICAL METHODS 



Here we describe how we numericalh 
We use mainly a relaxation technique % 



solve Eq. ©. 
7J, however, due 



to numerical problems at r = for g > 10 , we must 
combine this with a shooting method in order to achieve 
high numerical precision p8| . We begin by describing the 
relaxation method. 



A. Relaxation technique 

The idea of the relaxation method is to start from an 
ansatz, /, for a solution to the differential equation (DE), 
D(/) = 0, and then use a first order Taylor expansion to 
improve this ansatz iteratively. We want to find A/ so 
that D(/ + A/) = 0; expanding to first order we have 



D(/) + J dr 6 -^ • A/ = 



(27) 



Solving for A/ gives the improved function / + A/. To 
implement this on a computer, we first rewrite the ./V th 
order differential equation as N coupled first order equa- 
tions (in our case N — 2). We then put the system on 
a lattice and replace the differential equations by finite 
difference equations (FDE). The resulting algebraic equa- 
tions relate the values of the functions f = (f 1 , f 2 , f N ) 
(where f 1 = f, f 2 — df/dr...) at lattice points k (ffc) 
and k + 1 (ffc+i): 



D fe (r fe ,r fe+ i,ffe,f fe+ i} 







(28) 



here arc the FDEs at lattice point k. Taylor expand- 
ing as in (K7n gives 



D fc (f fc + Affc.ffc+i + Affc+i) « D fc (f fc) ffc+i) 



N 



jv 



(29) 



These equations are then solved for the increments Affc 
and the procedure is iterated until the corrections, Af^, 
are small enough. As a measure of the error we use: 



M N 



err 



M ■ N 



(30) 



k=l a=l 
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where M is the number of lattice points [M . If the error 
is large only a fraction of the correction is used when cal- 
culating the new function, this reduces the risk for over- 
correction near / = ±1 where the DE is very sensitive to 
disturbances in /. After each iteration we calculate the 
energy, using equation (||), and check that it decreases. 

B. Numerical solution for small g 

There are two problems to solve before we can apply 
the relaxation method: First, the boundary conditions 
are given at r = and at r = oo but we can only cover 
a finite interval [0, r max ] with the M lattice points (with 
the DE in the above form). This is a problem for large 
skyrmions, i.e. when g is small. Second, the relaxation 
method relies on the validity of a first order Taylor ex- 
pansion, therefore it is crucial to make an ansatz that is 
close to the true solution. 

To minimize the finite size effects we solve the equa- 
tions for several different r max and extrapolate to infinite 
system size. Fig. ^| shows how the energy scales with 
l/r mnl for g — 5 ■ 10~ 4 and 5 • 10 -6 . We see that the 
energy is more or less independent of system size at the 
largest sizes, whence the finite size effects are negligible. 
The spin, and the spin radius, show bigger finite size ef- 
fects and scaling to infinite system size is necessary to 
obtain accurate values when g < 10~ 6 . 
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FIG. 2. Skyrmion energies (per charge Z) as functions of 
system size for (a) g = 5 ■ 1CP 4 and (b) g = 5 • 1CF 6 , at v — 1. 

To obtain a good initial ansatz is a bit tricky; we pro- 
ceed as follows. The Belavin-Polyakov solution, Eq. (|l3|), 
to the scale invariant model is expected to be a good ap- 
proximation to the true solution in the core region (for 



some value of the scale parameter A), see Sec. As 
our ansatz we therefore use this solution with the func- 
tion const ■ r z added, where we choose the constant so 
that the boundary condition f(r max ) = 1 is fulfilled. We 
run the relaxation procedure for different A and choose 
the one for which the convergence is best. 

Having determined an ansatz we can start iterating. 
The Coulomb term Vc(r,f), which enters the DE, de- 
pends on the integral over /'. Thus we start each itera- 
tion by calculating the function Vc(r, /') using the most 
recent /. Then we iterate until the error is small enough, 
each time checking that the energy decreases. 

To check the accuracy of the solution we use a shoot- 
ing test. The relaxation solution is used as input to the 
shooting method (fourth order Runge-Kutta). We pick a 



G 



point ro close to r = (e.g. the first lattice point where 
/' 7^ 0) and treat the problem as an initial value problem 
with /(ro) and /'(ro) given from the relaxation solution. 
If the shooting solution thus obtained coincides with the 
relaxation solution we are confident that this is the right 
solution. Comparing the energy obtained from the relax- 
ation and the shooting solutions we find, for both Z = 1 
and Z = 2, that for the smallest g (g m 10~ 7 ) this dif- 
ference is of the order of 10~ 5 . The difference is approx- 
imately constant up to g = 10~ 4 — 10~ 3 . For g > 1CP 3 
the difference starts to grow; this is also reflected in the 
convergence properties of the relaxation method. This 
indicates that the ansatz we use deviates more from the 
true solution and the relaxation method has problems 
avoiding the singular behavior near r = 0. 
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C. Larger g 

The problems for larger g (g > 1CP 3 ) has to do with 
the structure of the differential equation ( [To| ) at the origin 
together with a smaller overlap between the solution and 
the ansatz. In this region we proceed as follows. The so- 
lution obtained from an initial set of relaxation iterations 
is used to shoot from the point ro p3] (close to r = 0) to 
a point a finite distance away from r = 0. This point is 
then used as a boundary point for a set of new relaxation 
iterations. With this method, which is illustrated in Fig. 
||, we avoid the problems near r = p0| ]. 

This method certainly introduces an error in the cal- 
culation. As a rough estimate of the error we take the 
energy difference between the initial relaxation solution 
and the solution obtained by the combined relaxation 
and shooting method. For g — 0.02 this is of the order 
of one percent. (Hartree-Fock theory Q shows that the 
skyrmion has a size of a few magnetic lengths and consists 
of a few flipped spins at g = 0.02. Whereas Hartree-Fock 
should be correct in this region, the validity of the ef- 
fective long-distance theory becomes questionable at dis- 
tances of the order of the magnetic length.) For a com- 
parison between Hartree-Fock and effective action results 
see Sec. below. 



FIG. 3. Sketch of the numerical method used to solve the 
equations of motion for g > 10 -3 . 



V. NUMERICAL RESULTS 



We have applied the technique described in Sec. [TV 
to skyrmions with topological charge Z = 1,2 and 3 at 
v = 1, i.e. v FM — 1, p s — l/16\/27r. Our results are 
presented in Figs. |||]. The Z — 3 skyrmion is not found 
to be the lowest energy excitation for any g considered 
here (in agreement with the analytic result, Sec. Ill), we 
therefore show results for Z = 1 and Z = 2 skyrmions 
only. 

Fig. U shows the z-component of the spin, n z = /, 
as a function of the radius, r. An analysis shows that 
for Z = 2 at g = 5 • 10~ 4 , the numerical solution agrees 
very well in the core region with the Belavin-Polyakov 
solution (|l3|), with A determined by Eq. (|l7j). For Z = 
1, no analytic solution is available for the values of g 
shown in the figure, however for g < g C3 = 4.3 • 10~ 5 , the 
numerical solution agrees well with the analytic solution 
(p"3|), with A determined by Eq. (p3|). Away from the 
core region, the relaxation solutions differ from the BP 
solutions, although the differences are quite small for Z = 
2. 
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FIG. 4. The relaxation solution, f(r), for g = 0.02, 0.0005 
and Z = 1, 2. Z — 2 is compared to the Belavin-Polyakov 
solution ([Til), with the scale, A, taken from Eq. (^). 

From the numerical solution f(r) we calculate the en- 
ergy of the skyrmion using Eq. (^). In Fig. ||, the energy 
of the Z = 1, 2 skyrmions, relative to the groundstate en- 
ergy £?(/ = i) = —jgr^axj ls given as a function of <7 and 
compared to Hartree-Fock (HF) results We see 

that the effective and Hartree-Fock theories agree very 
well for small g, but that the HF energy is substantially 
lower for larger values of g. In particular, the effective 
theory predicts that the skyrmions are the lowest en- 
ergy charged excitations for g < 0.018 whereas HF gives 
g < g Cl = 0.054. The difference between effective action 
and Hartree-Fock results is smaller for Z = 2 than for 
Z = 1. 
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FIG. 5. Effective action and Hartree-Fock results for the 
energy, E, per charge (relative to the ground state) of the 
Z = 1, 2 skyrmions as function of g. (Due to finite size effects 
in the HF calculation the HF energy becomes larger than the 
effective action energy for very small g.) The arrow in the 
inset marks the energy at g — 0. 

The effective theory is a long distance theory valid 
only for large smooth skyrmions, i.e. for small enough 
g. Hartree-Fock theory, on the other hand, is valid all 
the way from polarized quasiparticles and very small 
skyrmions near g Cl = 0.054 to quite large skyrmions 
(where finite size effects eventually limit the applicability 
of HF) . It is interesting to note that the effective theory 
works quite well also for fairly large values of g. In the 
limit of very small g, the effective theory becomes exact 
(see Sec. [II) and E/Z approaches yV/32 (marked by — > 
in the figure). The blowup in Fig. |^ shows a region (close 
to g = 0) where the Hartree-Fock energy is larger than 
the effective action energy, this is an effect of the diffi- 
culty of handling large enough systems in Hartree-Fock. 
Since in this range, the skyrmions become very large (see 
Fig. B) finite size effects make the Hartree-Fock ener- 



gies unreliable. For g < g c 



A x 10~ 5 , the Z 



skyrmions actually have lower energy (per charge) than 
the Z = 1 skyrmions, i.e. the Z = 1 skyrmions bind in 
pairs pa|. 



8 



We now compare the analytic small g expansions for 
the energy, derived in Sec. Ill, with the numerical effec- 
tive action results. Fig. ^ shows the skyrmion energies 
for 10~ 7 < g < 10~ 6 ; the inset shows the Z = 2 energies 
for larger g. For Z = 2 the agreement is excellent all the 
way up to g = 0.06. For Z = 1, the numerical results are 
seen to approach the analytic result Eq. (^4|) as g — > 0, 
however, deviations are still visible at g = 10~ 7 . The size 
of these is consistent with the size of the ignored correc- 
tion terms in Eq. (|2^). Note that the correction terms 
decrease extremely slowly with decreasing g. 
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FIG. 7. Spin, S z , for Z — 1,2 skyrmions; numerical ef- 
fective action results compared to HF and analytic small g 
expansions. The inset shows results for very small g where 
no HF data are available. For Z = 1, the analytic expansion 
breaks down at g C3 — 4.3 x 10~ 5 . 
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Fig. I] shows how the spin radii, R s , see Eq. (|2J 
of the skyrmions decrease with increasing g. The Z — 2 
spin radius is smaller than the Z = 1 radius for g < 10~ 3 . 



1*10 



FIG. 6. Numerical effective action results (points) for the 
energy compared with the analytic expressions, Eq. ([mJ) and 
([i"7|), (lines) for Z = 1, 2 skyrmions. The error bars come from 
finite size effects, for Z = 2 the error bars are smaller than 
the dots. The inset shows the Z — 2 data for large g. 

Fig. fj] shows the spin, S z , of the Z = 1,2 skyrmions 
relative to the groundstate, see Eq. (p5|). Since S z = 
this also gives the (topological) charge radius, R c . 
The numerical effective action results and the Hartree- 
Fock results become indistinguishable, for Z = 1 as well 
as for Z = 2, as g approaches zero but, in fact, agree 
quite well over the whole range although deviations are 
observable for larger g. The analytic small g expansions 
agree extremely well with the numerical effective action 
results for all g for Z = 2 and for g < g C3 for Z = 1 . 
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FIG. 8. Spin radii for Z = 1, 2 skyrmions; numerical ef- 
fective action results compared to HF and analytic small g 
expansions. The spin radii for Z = 1 and Z — 2 are equal at 
g~lCT 3 . 



In Fig. |9| we show spin and charge density pro- 
files. A comparison between the Hartree-Fock and ef- 
fective action densities reveals almost identical profiles 
at g = 5 • 10 -4 , see 32|, whereas at g = 0.02 the profiles 
deviate substantially. We also note that the agreement 
for the Z = 2 skyrmion is better than the agreement for 
the Z = 1 skyrmion (at large g). 
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FIG. 9. Spin and charge density profiles for Z — 1,2 
skyrmions; effective action and HF results. 



VI. ASYMPTOTIC SHAPE OF THE SKYRMION 

Here we will examine the detailed profile of the 
skyrmions in two regions, near the center and far from 
the center, where it is possible to compute it rigorously. 
This serves as a further check on our analysis. 

A. Small r 

Since the charge distribution is radially symmetric it 
is clear that dVc / dr\ r =o = 0. This implies that both the 
Zeeman and the Coulomb terms in the differential equa- 
tion ([h]) can be neglected and the solution is given by the 
Belavin-Polyakov solution, Eq. (|l3|). This is consistent 
with the numerical solutions, see Fig. |]. 



The crossover to this regime can be estimated by exam- 
ining the behavior of the Bessel functions for small values 
of their argument: 



K z (x) 



•>z-\ 



T{Z) 



(32) 



which shows that the crossover sets in at x z ~ 2 Z ~ 1 T(Z), 
for example, at r ~ R% (Z = 1) and r ~ V2R3 (Z = 2), 
i.e. farther out for higher Z. To compare these observa- 
tions with the numerical solution we plot ln[r(l — f(r))] 
against r for Z — 1 and Z = 2, at v = 1, see Fig. UXj . 
When g — 5 ■ 10~ 6 , R3 = 250 and we get a straight line 
for r > 400 when Z = 1 and for r > 700 when Z = 2; in 
qualitative agreement with our above observation. Read- 
ing off the slope for Z = 1 in the figure and comparing to 
Eq. (^lj) one finds R3 — 247 ± 9, in excellent agreement 
with the exact value R3 — 250. For Z = 2, the exponen- 
tial region sets in further out and, in addition, the curve 
bends up slightly for r > 1000 due to finite size effects in 
combination with the finite resolution of real numbers in 
the computer, therefore this curve is harder to analyze. 
However, comparing the curve for several different sys- 
tem sizes and keeping the part which is unchanged, one 
finds that it is consistent with an exponential behavior 
with R 3 = 250. 
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2000.0 



FIG. 10. Numerical effective action result for the asymp- 
totic behavior of f(r), at v = 1. Both the Z = 1 and the 
Z = 2 cases are well described by the asymptotic formula 
derived in the text. 



B. Large r VII. EDGES 



As we showed in Sec. far from the core, r > R 3 , i n previous work 0| we have used the effective theory 

the skyrmion profile, 1 - f(r) oc K%(r/R 3 ), is exponen- to study the reconstruction of a sharp polarized edge of a 
tiab semi-infinite quantum Hall ferromagnet as the confining 
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potential is softened. The problem of reconstruction is 
as follows. For hard confinement, the ferromagnetic bulk 
state continues all the way out to the edge and the den- 
sity drops abruptly to zero over the scale of the magnetic 
length. When the edge potential softens, charge even- 
tually moves out, i.e. the edge reconstructs. The ques- 
tion is how this happens given the competition between 
the correlation energy of the quantum Hall state and the 



Hartree energy due to the confinement. In 14 we argued 
that the initial instability, as the potential softens, is to 
forming a spin texture along the edge |L5| . Charge is then 
moved out as a consequence of the spin-charge relation 
of spin textures. Here we will give details, previously 
unpublished, of how the formation of edge textures is 
analyzed within the effective theory. (Before proceeding, 
the following caveat is in order. The use of an unmodified 
magnetic action near the edge is not unproblematic given 
that there are gapless density modes near the edge among 
other things. While we have not attempted to derive an 
action for the edge from microscopies, we have checked 
|l4| by comparison with Hartree- Fock calculations, that 
the action we use is accurate for getting static properties 
of the edge at small g.) 

We consider a semi-infinite Hall system occupying the 
region x < (y is the direction along the edge). As an 
ansatz for the spin vector field n(r) we take 



n-x = yl - P(x) cos(ky) , 
n v = \J\ - f 2 (x) sin(fcy) , 
n z = f(x) ■ 

This leads to the topological density 
q(x) = ±-f(x) , 



(33) 



(34) 



where f'(x) = df jdx. The z— component of the spin, 
f(x), gradually falls from one in the bulk, to some value 
f e (— 1 < f e < 1) at the edge. As one moves along the 
edge the spin in the xy-plane rotates with wave number 
k. f e and k are continuous parameters that are to be 
determined. Note that, although locally this texture is 
identical to the texture for bulk skyrmions there is no 
quantization: J d 2 rq(r) can take any value. 
The electron density, p(x), is 



p = p + v FM q(x) , -oo < x < x R 
p = , x R < x . 



(35) 



The unreconstructed, sharp polarized edge corresponds 
to taking f(x) = 1, x < and f(x) = 0, < x. This 
case is obtained by letting q(x) = 0, xr — in (35). 
When the edge is textured, charge is moved out and xr is 
determined by charge conservation: J*^ dxp — j ^ dxp. 
This gives, 



xr = T^ 1 ~ ^) 



(36) 



We assume the confining potential is caused by a pos- 
itively charged background density, pb{x), that falls lin- 
early from its bulk value, p — v FM /2Tr, to zero over a 
region of width w centered around x = 0. The strength 
of the confining potential is varied by varying the param- 
eter w. For small w the sharp polarized edge is favored, 
but as w increases the edge will reconstruct. 

For a spin texture of the form ( |33| ) one finds the fol- 
lowing gradient, Zeeman and Coulomb energy densities 
(per unit length of the edge) : 



fXR r,2 



'XL 
47T 



k r XR 

j(l-/ e 2 )-/ dxf 

4 J XL 



En = -Vi 



/2 r w/2 



dxdx 5 p(x)5 p{x') In \x — x'\ 



(37) 



Here, xl is chosen such that fix) = 1 for x < xl, and 
5p(x) = p{x) — Pb{x) is the deviation of the electron den- 
sity from the background density. 

Varying the total energy E — Eq + Ez + Eq with 
respect to /, we find 



f"(l-f)+ff 



[k 2 f- 



4irp s 



dV c /dx](l-f 



2\2 



(38) 



where 



rXR j i u 

Vc(x,f')= j —[ Pb (x')-p-—f(x')}\n\x~x'\ 



2tt 

" /2 dx> 



4tt" 



f w/ dx' 

/ —p b (x')\n\x-x'\ 
Jx R 2vr 



(39) 



We solve Eq. (B8h using the relaxation method de- 



scribed in Sec. IV. As an initial ansatz for /, we take a 
modified BP solution f BP (x) = {x 2 - 4X 2 )/(x 2 + AX 2 ) 
(c.f. Eq. ([l3"|)) that obeys the boundary conditions 
f(x L ) = 1, f(x R ) = U 



f ans (x) = f e +(l-f e ) 



1 + Ibp{xr - x) 
1 + f BP (xR - x h ) 



(40) 



The ansatz depends on one parameter A (for given f e ). 
In the numerical procedure we vary this parameter until 
we get convergence. 

We can now study whether a spin textured edge has 
lower energy than the sharp polarized edge for given pa- 
rameters (g,w). In particular, we can determine the 
phase transition line in the (g, w)— plane where the sharp 
edge reconstructs by forming a spin textured edge. Solv- 
ing Eq. d38| ) using the relaxation method for given 
boundary condition f e and wave number k determines 
f(x). The energy of this state is calculated from Eq. 
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( |37| ) and k is determined by minimizing this energy. This 
gives the energy of the state E(g, w, f e ). f e = l gives the 
energy for the sharp edge. We find that the spin textured 
edge is lower in energy than the sharp edge for any f e if 
and only if the derivative of the energy at f e = 1 is pos- 
itive: dE(g, w, /e)/<3/e|/ e =i > 0. The result is that the 
textured edge has lower energy in a triangular region of 



parameter space, for small g and big enough w, see |14 



VIII. SUMMARY 

We have used a long wavelength magnetic effective 
action to compute various properties of quantum Hall 
skyrmions and find that these compare favorably with 
those obtained by more microscopic Hartree-Fock calcu- 
lations and by analytic small g expansions. We have 
also used the effective action to reliably predict the onset 
of texturing at the edges of quantum Hall fcrromagnets 
for small values of the Zeeman energies. Taken together, 
these show that this formalism and the relaxational tech- 
nique presented here have great utility in the study of 
quantum Hall ferromagnets. 
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